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Abstract 

We present preliminary results of a search for B° — > K°K° decays in the K®K® final state using 
a sample of approximately 23 million BB pairs collected by the BABAR detector at the PEP-II 
asymmetric B Factory at SLAC. We find no evidence for a signal and set an upper limit on the 
branching fraction of 7.3 x 10~ 6 at the 90% confidence level. 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in 
the understanding of CP violation in the B system. Measurements of rates and CP-violating 
asymmetries for B decays into charmless two-body final states provide information on the angles 
a and 7 of the Unitarity Triangle. However, in contrast to the theoretically clean determination of 
the angle (5'va.B decays to charmonium final states Q , the extraction of CP- violation parameters in 
charmless decays is complicated by hadronic uncertainties that are difficult to calculate from first 
principles. Accurate branching fraction measurements provide critical tests of theoretical models 
that are needed to obtain reliable information on a and 7. 

The BABAR collaboration has recently published (2| measurements of the branching fractions 
for B meson decays to the charmless hadronic final states 7r + 7r~, K + ir~ , K + ir°, K°ir + and K°tt°, 
upper limits on the decays to 7r + 7r° and K°K + and the results of a search for charge asymmetries 
in the modes B° — > K + tt~, B + — > K + tt° and B + — > K°ir + ^\ In this paper we report preliminary 
results of a search for B° — > K°K° decays through detection of the K®K® final state. Although 
the decay rate for B° — > K°K° is expected to be small (10 -6 -10 -7 ) in the Standard Model Q, 
final state rescattering effects can lead to enhancement of the branching fraction and the possi- 
bility of large strong phases, with correspondingly large CP- violating charge asymmetries 0, Q. 
Such rescattering effects may also have consequences for constraints on 7 derived from B — > Kit 
decays ||. Observation of the K°K° decay mode would provide important information about the 
strength of final state rescattering in charmless B decays. 

2 Data Sample 

The data used in this analysis were collected with the BABAR detector at the PEP-II e + e~ stor- 
age ring. The sample corresponds to an integrated luminosity of 20.6 fb -1 taken near the T(4S) 
resonance ( "on-resonance" ) and 2.6 fb -1 taken at a center-of-mass (CM) energy 40MeV below the 
T(45*) resonance ("off-resonance"), which are used for continuum background studies. The on- 
resonance sample corresponds to (22.57 ± 0.36) x 10 6 BB pairs. The collider is operated with 
asymmetric beam energies, producing a boost (^7 = 0.56) of the T(4S) along the collision axis. 
The boost increases the momentum range of two-body B decay products from a narrow distribution 
centered near 2.6GeV/c to a broad distribution extending from 1.7 to 4.3GeV/c. 

The BABAR detector is described in detail in Ref. ||. The primary detector element used in 
this analysis is the tracking system, which consists of a 5-layer, double-sided, silicon vertex detector 
and a 40-layer drift chamber filled with a gas mixture of helium (80%) and isobutane (20%). Both 
tracking detectors operate within a 1.5 T superconducting solenoidal magnet. 

3 Event Selection and K® Reconstruction 

Hadronic events are selected based on track multiplicity and event topology. Backgrounds from 
non-hadronic events are reduced by requiring the ratio of Fox- Wolfram moments H2 / Hq [j?]] to be 
less than 0.95 and the sphericity [|| of the event to be greater than 0.01. 

Candidate K® mesons are reconstructed in the tt + it~ final state from pairs of oppositely charged 
tracks that form a well-measured vertex and have an invariant mass within 11.2MeV/c 2 (3.5o") of 

1 Charge conjugate modes are assumed throught this paper. 
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Figure 1: Invariant mass of 7r + 7r~ pairs in the inclusive Kg sample selected with requirements on 
the angle between the flight direction and momentum of the K®, and on the transverse momenta 
of the decay products relative to the K® momentum. 

the nominal K® mass ||. The measured proper decay time of the K® candidate is required to 
exceed 5 times its error. Figure |l| shows the invariant mass distribution of an inclusive sample of 
high momentum (> lGeV/c) K® candidates. 

4 B Reconstruction and Background Rejection 

B meson candidates are reconstructed by combining pairs of K® candidates. The kinematic con- 
straints provided by the T(4S) initial state and relatively precise knowledge of the beam ener- 
gies are exploited to efficiently identify B candidates. We define a beam-energy substituted mass 
m-Es = \J — p 2 B , where E\> = (s/2 + p« • ps)/^, yfs and Ei are the total energies of the e + e~ 
system in the CM and lab frames, respectively, and pj and pb are the momentum vectors in 
the lab frame of the e + e~ system and the B candidate, respectively. The mEs resolution for B 
decays into all-charged final states is dominated by the beam energy spread and is determined 
to be 2.6 ± 0.1MeV/c 2 from a Gaussian fit to a large sample of fully reconstructed B decays. 
Candidates are selected in the range 5.2 < ttt-es < 5.3GeV/c 2 . A sideband region is defined as 
5.2 < m ES < 5.26GeV/c 2 . 

We define an additional kinematic parameter AE as the difference between the energy of the 
B candidate and half the energy of the e + e~ system, computed in the CM system. The AE 
distribution for signal events is a Gaussian centered near zero. The resolution on AE is estimated 
to be 21 ± 5 MeV based on Monte Carlo simulated B° — > K®K® decays and the observed difference 
in widths between a control sample of fully reconstructed B decays in data and in Monte Carlo 
simulation. We require |A.E| < 0.1 GeV. A sideband region is defined as 0.1 < \ AE\ < 0.3 GeV. 

Detailed Monte Carlo simulation, off-resonance data, and events in on-resonance sideband re- 



40 I 




-4 -2 2 4 -4 -2 2 4 

Fisher variable Fisher variable 



Figure 2: Distributions of Fisher discriminant output for Monte Carlo simulated B — > K^Kg 
decays (left) and events in the on-resonance toes sideband data (right). 

gions are used to study backgrounds. The contribution due to other B meson decays, both from 
b —* c and charmless decays, is found to be negligible from a detailed Monte Carlo study. The 
dominant source of background is from random combinations of true K® mesons produced in the 
e + e~ — > qq continuum events (where q = u, d, s, or c). In the CM frame this background typically 
exhibits a two-jet structure, in contrast to the isotropic decay of BB pairs produced in T(4S) 
decays. We exploit the topology difference between signal and background by making use of two 
event-shape quantities. 

The first variable is the angle 9 S between the sphericity axes of the B candidate and of the 
remaining tracks and photons in the event. The distribution of | cos 9 S \ in the CM frame is strongly 
peaked near 1 for continuum events and is approximately uniform for BB events. We require 
| cos# s | < 0.9, which rejects 66% of the background remaining at this stage of the analysis. 

The second quantity is a Fisher discriminant T constructed from the scalar sum of the CM 
momenta of all tracks and photons (excluding the B candidate decay products) flowing into nine 
concentric cones centered on the thrust axis of the B candidate. Each cone subtends an angle of 
10° and is folded to combine the forward and backward intervals. Monte Carlo samples are used 
to obtain the values of the Fisher coefficients, which are determined by maximizing the statistical 
separation between signal and background events. Figure § shows distributions of T for Monte 
Carlo simulated B° — ► K®K® decays and background events in the toes sideband region. No 
cut is applyed on J 7 , instead the distributions for signal and background events are included in a 
maximum likelihood as described in the next section. 

A total of 286 candidates in the Runl on-resonance data satisfy our selection criteria ( | cos 9 S \ < 
0.9, 5.2 < toes < 5.3GeV/c 2 and \AE\ < 0.1 GeV) and enter the maximum likelihood fit. The 
total selection efficiency is (36.6 ± 4.6)%, where the error is dominated by uncertainty on the K® 
reconstruction efficiency (6% relative error per K®). 



5 Signal Extraction 



Signal and background yields are determined from an unbinned maximum likelihood fit using m^s, 
AE and T. The likelihood is defined as 

N 

£ = e ~(N s+ N B ) "Q (n sV S (mES) AE 1 , J 7 ) + iV B pf (m ES , AE..F)) , (f) 

i 

where N$ and are the fitted number of signal and background events, respectively; iV is the total 
number of events entering the fit; and V{ is the product of probability density functions (PDFs) 
for mES) AE, and T that are assumed to be uncorrelated. The quantity — 21n£ is minimized with 
respect to the fit parameters N$ and Nb- 

The parameters for background ?ties and AE PDFs are determined from events in on-resonance 
AE and /ties sideband regions, respectively. The m-ES shape is parameterized by a threshold 
function (l(J /(mEs) w-esv 7 ! — x 2 exp[— £(1 — x 2 )], where x = mEs/ m o and m o is the average 
CM beam energy. The background shape in AE is parameterized as a second-order polynomial. 
The signal distributions for mES an d AE are described by Gaussians, where the mES mean and 
width are determined from a sample of fully reconstructed B decays while the AE mean and 
width are estimated from Monte Carlo simulated B° — > K®K® decays and scaled according to 
the observed difference between a control sample of fully reconstructed B decays in data and in 
Monte Carlo simulation. Events in on-resonance rriEs sideband regions and Monte Carlo simulated 
signal decays are used to parameterize as the sum of two Gaussians the Fisher discriminant PDFs 
for background and signal. Alternative parameterizations for T , obtained from off-resonance data 
(for background) and fully reconstructed B decays (for signal), are used to determine systematic 
uncertainties. 

The fitted number of signal events is Ng = 3.4lt|; 4 , where the error is statistical only. Figure |3| 
shows the tties and AE distributions for events satisfying the selection criteria and additional 
requirements on the likelihood ratio NgVs/iNs'Ps + NbVb), where the probabilities include all 
variables (wies, AE, J 7 ) except the one being plotted. The likelihood ratio cuts are chosen to 
minimize the upper limit on the branching fration, and the curves represent the fit result scaled by 
the efficiency of the additional requirements. The statistical significance of the signal yield is 1.5, 
calculated as the square root of the change in — 21n£ when the yield is fixed to zero. We conclude 
there is no evidence for a signal and calculate a 90% confidence level Bayesian upper limit as the 
value of Atjl for which J^ UL C max dNs / Anax dN$ = 0.90, where £ max is the likelihood as a 
function of N$, maximized with respect to the remaining fit parameter (Nb)- The resulting upper 
limit on the yield is A%l = 9. The fitting procedure has been validated with extensive Monte Carlo 
studies. 

6 Branching Fraction Results 

The branching fraction is defined as 

e is the total KgK® selection efficiency, N BB - = (22.57 ±0.36) x 10 6 is the total number of BB pairs 
in the dataset, and B(K® — * ir + ir~) = 0.6861 ||. We assume the Standard Model prediction that 



> 

<u 



^ 6 



a 

u 
> 



BaBar 



> 

o 



fl 3 

> 



BaBar 



2 - 



5.2 



5.225 



5.25 
m ES 



5.275 



5.3 



(GeV/c ) 



-0.1 



-0.05 





AE 



0.05 0.1 

(GeV) 



Figure 3: Distributions of mES (left) and AE (right) after additional requirements on likelihood 
ratios. The curves represent projections of the maximum likelihood fit result. 

B° — > K®K® proceeds through the K°K° intermediate state (as opposed to K°K° or K°K°) and 
use B(K°K° — » K®K®) = 0.5.0 Implicit in Eq. |2| is the assumption of equal branching fractions for 
r(45) -> and r(4S) -» 5+B~ 

Systematic uncertainties on the branching fraction arise primarily from uncertainty on the 
reconstruction efficiency and uncertainty on N$ due to imperfect knowledge of the PDF shapes. 
The latter is estimated either by varying the PDF parameters within la of their measured uncer- 
tainties or by substituting alternative PDFs from independent control samples. Table [l] summarizes 
the various sources of systematic error on N$, where the total error is calculated as the sum in 
quadrature of the individual uncertainties. 

We measure a central value branching fraction of B(B° -»■ K°K°) = (1.8t\\t ± 1-8) X 1(T 6 , 
where the first error is statistical and the second is systematic. An upper limit is calculated 
by increasing N\jh and decreasing the efficiency by their respective systematic errors. We find 
B(B° -► K°K°) < 7.3 x 10" 6 at the 90% confidence level. This result is a significant improvement 



over the existing upper limit from the CLEO Collaboration [11|, and is approaching the upper 
range of current theoretical estimates. 
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Table 1: Summary of absolute systematic errors (<7jv s ) on the signal yield due to imperfect knowl- 
edge of the PDF shapes. The total error is calculated as the sum in quadrature of the individual 
uncertainties. 
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